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We argue that from observations alone, only the transverse power spectrum Ci{z\,zv) and the 
corresponding correlation function z\, Z2) can be measured and that these contain the full three 
dimensional information. We determine the two point galaxy correlation function at linear order 
in perturbation theory. Redshift space distortions are taken into account for arbitrary angular 
and redshift separations. We discuss the shape of the longitudinal and the transversal correlation 
functions which are very different from each other and from their real space counterpart. We then 
go on and suggest how to measure both, the Hubble parameter, H(z), and the angular diameter 
distance, Da(z), separately from these correlation functions and perform an Alcock-Paczynski test. 

PACS numbers: 98.80.-k,98.80.Es 
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I. INTRODUCTION 

Cosmology has become a data driven science. After 
the amazing success story of the cosmic microwave back- 
ground (CMB), see [TH3], which is still ongoing [4], we 
now also want to profit in an optimal way from actual 
and future galaxy catalogs. Contrary to the CMB which 
comes from the two dimensional surface of last scattering, 
galaxy catalogs are three dimensional and therefore con- 
tain potentially more, richer information. On the other 
hand, galaxy formation is a complicated non-linear pro- 
cess, and it is not clear how much cosmological informa- 
tion about the underlying matter distribution and about 
gravitational clustering can be gained by observations of 
the galaxy distribution. This is the problem of biasing 
which we do not address in this paper. Here, we sim- 
ply assume that on large enough scales biasing is linear 
and local, a hypothesis which might turn out to be too 
simple [5]. 

When observing galaxies we measure their redshift 
and angular position. To convert this into a three- 
dimensional galaxy catalog we must make an assumption 
to relate the observed redshift to a distance. For small 
redshift, the simple relation D = z/Hq can be used. Red- 
shift space distortions (RSD) can be taken into account 
with a convenient expansion in tripolar spherical har- 
monics O [7] . This gives an accurate description of the 
correlation function at small scales. Apart from RSD, 
a wrong measurements of Ho will just rescale the entire 
catalog but not distort its clustering properties. 

However, if wc go out to high rcdslrifts, z > 1, non- 
linear terms in z become relevant, and wrong assump- 
tions about the distance redshift relation can bias the 
entire catalog. We therefore believe that it is important 
to analyze the truly observed catalog, either using the 
Ci(z\,Z2) spectra introduced in Ref. [8 , or the angular 
correlations functions Zi,z 2 ) to describe the observa- 
tions, and to compare them with their theoretically ob- 
tained counterparts. In this way we truly compare obser- 
vations with their theoretical modeling. If, on the other 
hand, we determine a power spectrum in Fourier space 



for the observed catalog, P(k), we have already assumed 
a cosmology to convert observable redshifts into length 
scales. Therefore, e.g., cosmological parameter estima- 
tions using P(k) can at best be viewed as a consistency 
check. If the cosmological parameters used to obtain 
P(k) agree with those revealed in a Markov Chain-Monte 
Carlo parameter estimation, the model is consistent. A 
thorough error estimation seems however, quite tricky. 

We therefore advocate to abandon this 'mixed' method 
in future, high redshift catalogs in favor of the more direct 
procedure which compares theoretical models with the 
directly observed two-point statistics, Ci(zi,z%) and/or 
£(6,z 1 ,z 2 ). 

This paper is structured as follows: in the next sec- 
tion we explain how to compute £(#,£1,22) an( A its co- 
variance matrix from the theoretical power spectrum. In 



Section III we discuss the longitudinal and transverse cor- 
relation functions and we show how the baryon acoustic 
oscillations in these correlation functions can be used as 



an Alcock-Paczyhski test (5]. In Section IV we conclude. 
In Appendix [A] we generalize the expansion in tripolar 
spherical harmonics of [7] to model the RSD of the cor- 
relation function at arbitrary redshifts. 



II. GENERALITIES 

In an observation which simply counts galaxies, we 
measure the galaxy distribution in angular and redshift 
space. For simplicity, and since we are interested only in 
large scales, we assume that its fluctuations are related by 
a scale independent bias factor b(z) to the matter density 
fluctuations 



A(n,z) = 



p(n, z) - p b (z) 
Pb(z) 



Here p(n, z) is the true matter density in direction n from 
the observer which we position at x = 0, at measured 
redshift z and Pb(z) is the matter density of a Fricdmann 
background Universe. 
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In Ref. [5], A(n, z) has been computed including all rel- 
ativistic effects in first order in perturbation theory. In 
the present paper we include only the terms which con- 
tribute significantly (more than 1%) under the circum- 
stances discussed here. These are the density fluctuation, 
the redshift space distortion and, in some cases also, the 
lensing term. For simplicity, we shall not consider lens- 
ing in this work, see however [TUj. Neglecting the other 
gravitational and Doppler terms, the expression given in 
Ref. [8] for A(n, z) becomes 1 

A(n,z) = D(n,z)-~(V(r(z)n,t(z))-n) 

-^V(r(z)n,t(*)).n. (1) 

Here t(z) is conformal time, D is the density fluctuation 
in comoving gauge, V is the peculiar velocity in longitu- 
dinal gauge, and H = aH is the comoving Hubble param- 
eter. Setting a(z) = 2/H corresponds to the sub-leading 
term of the redhsift space distortion, but one can also 
consider other terms than RSD in the function a(z), like 
the boost — 'H/(Y'H 2 )V • n calculated in [S], or a selec- 
tion function [TT] . For more details see and [5] . 
The power spectrum P(k) of density fluctuations (e.g. 
in synchronous gauge) can be calculated with standard 
Boltzmann solvers [T31 Q3] . 

We define the redshift dependent angular correlation 
function by 

£(9,zi,z 2 ) = (A(n, 2i)A(n',z 2 )) with n n' = cos9 . 

(2) 

Statistical isotropy implies that £ depends only on the 
angle 9 and not on the directions n and n' separately. In 
Ref. [5] , this correlation function is expanded in spherical 
harmonics, 

t(9, zi, z 2 ) = i J2( 2£ + ^)Ci(z 1 ,z 2 )P i (cos(9)) (3) 
i 

where Pi(p) are the Legendre polynomials and we call 
Gf,(z\,z 2 ) the angular power spectrum. In principle one 
can now go on and expand the z\ and z 2 dependence in 
terms of orthonormal functions. This direction has been 
explored in Ref. |15j . Since the background Universe and 
therefore also the correlations depend on Z\ and z 2 sepa- 
rately and not simply on r(zi)—r(z 2 ), we refrain from this 
step here. Clearly, £{9,Z\,z 2 ) or the Cg(zi, z 2 )'s contain 
the full three dimensional 2-point clustering information, 
exactly like the power spectrum. However, the distinct 
advantage of these quantities w.r.t the commonly used 
power spectrum is that we need not assume any back- 
ground cosmology in order to infer them from the obser- 
vations. 



1 Note that in [8] n denotes the direction of propagation, which is 
opposite to the direction of observation considered here. 



In our approximation A(n, z) contains two terms which 
contribute to £ with their correlators. We indicate them 
as D for the density term and z for the redshift space 
distortion, 

£ sal (e,zi,z 2 ) = &i6a£i>2>(Mi>*&)+&i&J*(Mi,*a) 
+b 2 £ Z D(0, zt,z 2 ) + £ zz (9, z u z 2 ), (4) 

where hi = b(zi) is the bias which only multiplies the 
density but not the velocity. We assume negligible ve- 
locity bias. We now compute these terms within lin- 
ear cosmological perturbation theory. For this we de- 
note by P(k) the density-density power spectrum today 
and by G(z) the growth rate, such that G(z = 0) = 1 
and P(k,z) = G 2 (z)P(k). This product ansatz is justi- 
fied under the assumption of vanishing sound speed (pure 
matter fluctuations) such that all modes evolve with the 
same growth function. It is certainly sufficient in stan- 
dard ACDM at z ~ 10, but might have to be revised 
for different dark energy models or if massive neutrinos 
are taken into account [16] . If the neutrino masses are of 
the order of the measured mass differences, they are not 
very relevant for structure formation, but if their mass is 
larger, they lead to damping of small scale power. More 
precisely P(k) and G[z) are defined by 

(D(k,t )D*(k',t )) = (2^) 3 <S(k-k')P(fc) (5) 
D(k,t(z)) = G(z)D(k,t Q ). (6) 

The comoving distance r(z) from the observer to redshift 
z in a fixed cosmology is given by 

r dz ' 

This distance is related to the angular diameter distance 
via Da(z) = r(z)/(l + z) and to the luminosity distance 
via D L {z) = r(z)(l + z). 

Assuming vanishing spatial curvature, K — 0, the law 
of cosines implies that the comoving distance r(z%, z 2 , 9) 
between two points at arbitrary redshifts z\ and z 2 seen 
under an angle 9 is 

r(z u z 2 , 9) = y/r( Zl ) 2 + r(z 2 ) 2 - 2r(z 1 )r(.2 2 ) cos 9 . (8) 

With this we can express the angular correlation function 
£(9, zi, z 2 ) in terms of the spatial correlation function 
£(r, zi,z 2 ), 

£(9,z 1 ,z 2 ) = £(r(zi,z 2 ,9),zi,z 2 ) . (9) 

This is useful since £ is simply the Fourier transform of 
the power spectrum. For ^dd we now obtain 

S, DD (9,z ll z 2 ) = G l G 2 J — fc 2 P(fc)j o (fcr(z 1 ,z 2 ,0)), 

(10) 

where Gj — G(zj). In a real galaxy survey, we observe 
a number of galaxies within a certain redshift bin and 
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within a solid angle. We assume that on sufficiently 
large scales, the galaxy number fluctuations are related to 
the underlying density fluctuation by a scale- independent 
bias factor b(z) and the velocities are unbiased. Using the 
continuity equation, the redshift space distortion term 
can also be expressed in terms of the power spectrum P 
and we find 

C g ai(ni,n 2 ,z 1 ,z 2 ) = (A ga i(ni,Zi)A* al (n 2 ,2;2)) = 
d 3 k 



bidb 2 G 2 
1 



(2tt) 3 
Pi , 2ft 



P(k)e 



k r 



P 2 (n 1 -k) — Pi(ni-k) 

3 3 T\k 

/3 2 2/3 2 t ia 2 H 2 (3 2 „ , 

1 + y + ^- P2(n2 - k) + ^ 2 ^ Pl(n2 ' k) 



(11) 



where k = k/fc and Pj{x) is the Legendre polynomial 
of degree j. We also use r = ri — r 2 , oti = a(zi), 
Hi = U(zi), and fa = f{z i )/b{z l ), where f(z) = 
d log G(a(z)) / d log a(z) , a(z) = 1/(1 + z) being the scale 
factor. The second term in the brackets of equation ( 11 ) 



contributes to the well known plane-parallel approxima- 
tion of RSD [TTJ HH] . From now on we neglect contribu- 
tions to the function a(z) other than the RSD-rclatcd 
term 2/H, so that the third and fourth terms in the 
brackets just correct the plane-parallel approximation 
and render the formula valid for arbitrary separation an- 
gles. 

We define the functions 



dk 
2^2 



k m u{kr)P{k) , 



(12) 



where ji{x) denotes the spherical Bessel function of index 
£. The results are obtained in terms of such integrals with 
< t < 4 and < m < 2 multiplied by functions of the 
redshift and of 9. This can be expressed in (relatively) 
compact form using an expansion in tripolar spherical 
harmonics [6l [7] . We present the details in Appendix [A] 
Here we just explain the basic idea: it is useful to in- 
troduce coordinates where the triangle consisting of the 
observer and the two positions given by S± = (ni, z\) and 
5 2 = (n 2 ,z 2 ) lies in the plane d = ir/2 (the x-y plane) 



and the direction (S^Si) is parallel to the x-direction. 
We also define Si as the position with the smaller az- 
imuthal angle, hence <f>i < 2 . This is depicted in Fig.[T] 

Then 6 = fa — <f>i, and the angles fa, fa can be ex- 
pressed in terms of the observables z\ , z 2 and 9 by 



— sin t 
r 



9 2 — 9 = sin 



'2 



(13) 



Here Tj = r(zj) and r is the distance given in Eq. ([8]). The 
constraint fa < fa can be inverted using the symmetry 




FIG. 1: Coordinate system. The triangle formed by the ob- 
server O and the two galaxies is orthogonal to the direction 
of observation z. We assume fa < fa. 



£gal(n 2 ,ni,n,r) = ^g al (-n x ,-n 2 ,n,r). Using this coor- 
dinate system, the tripolar expansion presented in the 
Appendix only contains a small number of terms and we 
end up with an expression of the form 



£ ga l(zi> . ^2, 9) = b(zi)G(zi)b(z 2 )G(z 2 ) 

> [a ni n 2 cos(nifa) cos(n 2 fa) 



E 

ni, "2=0,1,2 



+b ni7l2 sm(nifa) sm(n 2 fa)} . (14) 



In a true observation we cannot measure the number of 
galaxies at a precise redshift, but we actually work with 
redshift bins around the fiducial redshifts Z\, z 2 , where 
the probability distribution around the redshift z is given 
by a window function, typically a Gaussian 



W(z, z) oc exp 



1 fz- z 
"2 \ W 



(15) 



with standard deviation <r z , centered around the mean 
redshift z and normalized to 1. The observable angular 
correlation function is then given by 



dz\W(z\, zi)b{zi)G{zi) x 

dz 2 W(z 2 ,z 2 )b(z 2 )G(z 2 ) x 

E(fa(zi,z 2 ,9),fa(zi,z 2 , 9),r{zi,z 2 , 9)) , 

(16) 



where we have introduced 

fa, r) = 



[a ni n 2 cos(ni^i) cos{n 2 fa) 

ni, 712—0,1,2 

+b niTl2 sin(mfa) sin(n 2 fa)} . (17) 
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The difference between £ ga i and £ s is that the latter has 
been smeared over redshifts with a window function of 
widths a z . The non- vanishing coefficients a niiTl2 and 
b ni n 2 are given in Eqs. (A7| and (A8| in terms of the 



and (i(x). With this, eq. ( 16 ) demands the same numer- 



integrals Q a . 

In order to understand the signal in clustering anal- 
ysis of galaxy surveys, we must also estimate the error 
bars in the measurements. Since we are interested on 
scales where non-linear effects do not play a prominent 
role, we can model the expected errors following [!§] . 
The covariance in the measurements of £ s (#, Zi, z 2 ), i.e., 
Cov w = (f s (6>, z u z 2 )£ s {0', z\,z%)), where | g denotes the 
estimator used for £ g , can be related to the angular power 
spectrum, defined in equation ([3]), as 

Cov cc , = ^^^±Ip,(cos0)P,(cos0')(^ s + l/7V) 2 

= ^^^PKcos0)P,(cos0')Cov„(z,z'). (18) 

This includes the effects of sampling variance, shot-noise 
(N is the number of objects per steradian), partial sky 
coverage (f s ky is the observed fraction of the sky), pho- 
tometric redshift uncertainties and RSD. We have intro- 
duced the power spectra smeared over a redshift bin. 

C?(zi,z 2 ) = J dz 1 W(z 1 ,zi)b{zi) x 

dz 2 W(z 2 ,z 2 )b(z 2 )C e (z 1 ,z 2 ) (19) 



This angular power spectrum has been discussed in [5] 
and, according to their notation, we are approximating 
C e = CP D +Cf z +Cl z , where the three terms refer to the 
correlations between density-density, density- RSD, and 
RSD-RSD, respectively. We emphasize that if we write 
the covariance matrix in terms of reconstructed distances 
instead of observable angles and redshifts, we would not 
obtain a consistent error estimation. Also, the covariance 
matrix of an ideal experiment is diagonal in I but not in 
z and z' . This complicates the analysis as compared, e.g. 
to CMB spectra. 

We shall consider the detailed computation of this co- 
variance matrix for a given experiment in future work, 
whereas here we focus on the correlation function. 



From a numerical point of view, equation ( 16 1 may look 



quite cumbersome, since it involves a sum of three dimen- 
sional integrals, two of which over redshift and one over 
wavelengths, see eq. (12 1. However, the functions (P(x) 



depend only on the cosmology so, given the cosmological 
parameters, they can be calculated once and then stored. 
Therefore, eq. ( |16[ ) requires only a two-dimensional inte- 
gral of the window functions over redshifts, which can be 
performed rapidly. Then, eqs. (A7), which are the lead- 
ing coefficients at BAO scales (see Fig. [3] below), require 
the evaluation of the only three functions (q(x) (which 
is also called the 'real space' correlation function), ( 2 (x) 



ical effort as the plane-parallel case, where the evaluation 
of three functions similar to £™(ir)'s is needed [18 . How- 
ever, eq. ( [To] ) takes fully into account wide-angle effects 
that are neglected in the plane-parallel approximation. 



To evaluate eq. (16), we compute the power spectrum 
using the Camb code [13] • The linear result is used to 
obtain the £|™(r)'s, equation (12), on large scales r > 



150 Mpc/h, while Halofit 20J~is used on smaller scales 
to take into account corrections for non-linear evolution 
approximately, as in |21) . 

In the examples presented in this paper we use the 
cosmological parameters TO = 0.24, ilf, = 0.04, 57a = 
0.76, h — 0.73, primordial amplitude and spectral index 
equal to A s = 1.99 • 10~ 9 , n s = 0.96, respectively, and a 
constant bias 6 = 1. 

r(z,z,8) [Mpc/h] 
100 200 300 400 500 600 700 



0.004 



0.002 



0.000 



-0.002 




FIG. 2: The angular correlation function |16f at z — 0.7 mul- 
tiplied with 6 2 (in degrees) is shown. The different curves cor- 
respond, from top to bottom, to the result in redshift space 
discussed above (top solid line), the plane- parallel approxi- 
mation (middle dashed line) and the real space correlation 
function (bottom, short-dashed line), respectively. 

In Figure [2] the angular correlation function integrated 
over redshift bins, Eq. (16), is shown. A Gaussian se- 



lection function compatible with photometric surveys 
a z /(l + z) = 0.03 is employed, and the bins are cen- 
tered around the mean redshifts z± = z 2 = 0.7. The top 
axis has been evaluated using equation Q. For better 
visibility of the BAO peak we plot 6» 2 £ s (0). The results 
in redshift space differ significantly from the real space 
correlation function. The plane-parallel approximation 
[H equation (10)], does not give an accurate description 
for angles 8 > 1°, compared to the wide-angle RSD result 
obtained with equation (14). 



Another interesting observation is that once we cor- 
rectly include redshift space distortions, the correlation 
function no longer passes through zero. The enhanced 
clustering observed in redshift space, has a correlation 
function which is everywhere positive at fixed observed 
redshift. For the pure real space correlation function this 
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behavior is not possible due to the integral constraint, 
r 2 i DD {r)dr = 0. 



III. 



THE ALCOCK-PACZYNSKI TEST 



As the correlation function is no longer isotropic when 
accounting for RSD, this constraint no longer holds for 
the observed £ s . The large-scale anti-correlation in real 
space is due to the initial conditions set up at the end 
of inflation and evolution under gravity. The transition 
of a Harrison-Zeldovich from P(k) oc k on large scales to 
P{k) oc k~ 3 spectrum is due to the fact that density per- 
turbations can grow only in the matter era. For modes 
that enter the horizon during the radiation dominated 
era (see, e.g., [3]), this delays the growth roughly un- 
til matter-radiation equality, so that k 3 P(k) is roughly 
constant on these scales. The scale of the horizon at 
matter-radiation equality is imprinted as the location of 
this turnover in the power spectrum [22] and equivalently, 
as the location of the zero-crossing of the real space cor- 
relation function [23] , 

However, along the transverse direction, RSD suppress 
anti-correlations at large scales. In fact RSD introduce 
additional positive correlations between velocity and the 
density contrast at fixed observed redshift. These addi- 
tional positive correlations are larger than the negative 
density correlations. This is shown in Figure [2j where 
a relatively large anti-correlation at 6? m i n ~ 6° in the 
real space correlation function £ DD , |^Q D (6' m i n ) #„ lin | ~ 
2.5 x deg 2 , is lifted to a small positive correlation in 
redshift-space |£ s (6Un) CiJ ~ 7 - 5 x 10~ 4 deg 2 . 
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FIG. 3: Residuals of the redshift space correlation functions 
calculated as in Figure [2] neglecting (i.e., considering only 
leading RSD contribution) and taking into account the sub- 
leading coefficients in equation (A8l for a(z) = 2/H(z). 



The curve in Figure [3] suggests that the full RSD result, 
with a(z) = 2/W(z), differs by up to 6% from the leading 
contribution (i.e., with a(z) = 0) only for large angles. 
However, on angular scales up to 9 ~ 7°, which is well 
above the BAO's peak, the difference is less than 1%. 
This justifies to neglect a(z) in the next section, where 
we are interested in the BAO scale. 



One of the deepest enigmas of modern cosmology is 
the problem of dark energy: The energy density in the 
present Universe seems to be dominated by a cosmolog- 
ical constant A providing about 70% to the expansion 
rate of the Universe. For a review of the dark energy 
problem see, e.g. [23]- So far, most indications for dark 
energy come from measurements of the distance redshift 
relation |25j . and therefore rely on the Friedmann equa- 
tion and on relation This relation can be tested if 
we can measure both, r(z) and H(z) independently. Al- 
cock and Paczyhski have proposed such a measurement 
as follows: imagine a spherical object of comoving size 
L in the sky at redshift z. The redshift difference of its 
front and back is then given by Az^(z) = LH(z) and 
we see it under an angle 0l{z) = (i+ z )d a (z) • Knowing 
L we can in principle determine both (1 + z)Da(z) and 
H(z) by measuring Azl and 0l- But even without any 
knowledge of L we can infer the product 

F(z) = (1 + z)H(z)D A (z) = = F AP {z) . (20) 

Combining this with a measurement of the luminosity 
distance Dl(z) = (1 + z) 2 Da(z), e.g. from supernova 
type la data [35], we can break the degeneracy between 
H(z) and Da(z). This allows us to test the relation 



D A (z) 



dz' 



z + lj H(z') 



or F(z) 



j H{z) 



dz' 
o H(z') 



which must be valid, if the geometry of our Universe is 
close to a Friedmann-Lemaitre metric. 

Assuming an incorrect cosmology for the determina- 
tion of the galaxy power spectrum causes geometric 
redshift-distortions, in addition to the dynamical redshift 
distortions due to peculiar velocities of the galaxies dis- 
cussed here. Hence, the observed clustering signal can 
be used to constrain the underlying cosmology (Alcock- 
Paczyhski test [9]). This is also why it is important to 
understand the dynamical RSD. As suggested also, e.g., 
in [27] . we want to propose the baryon acoustic oscilla- 
tions (BAO) scale tbao hi the galaxy correlation function 
to provide the scale L. 

When considering physical distances, a cosmology 
must be assumed and, e.g., iterative methods have been 
proposed to converge to the correct cosmology j2HJ- The 
inferred galaxy clustering in a different cosmological 
model can also be obtained from the fiducial one by a 
rescaling of the transverse and parallel separations [29j . 
However here, instead of using these approximate proce- 
dures, we shall not assume a cosmology to determine the 
power spectrum but we work only with the directly ob- 
servable redshift dependent angular correlation function 

In this work we discuss how to constrain the quan- 
tity F(z) = (1 + z)H(z)Da(z). The degeneracy be- 
tween Da{z) and H(z) can then be lifted e.g. when com- 
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paring this with the direction averaged BAO scale [50] . 
D\(z)/H(z), or with SN la distances. 

For an introduction to BAO's see [3TJ [32]. It has 
been suggested before to use the BAO's for an Alcock- 
Paczyhski test, e.g., in (33J [33]. The BAO ring (ob- 
served in [35]) has been studied in [33 [37], and [35] 
analyses the anisotropics affecting the BAO's without as- 
suming a particular form of RSD. Experimental realiza- 
tions of the Alcock-Paczyhski test analysing the cluster- 
ing anisotropies are reported in 39 -|T]. In [37J the pos- 
sibility to observe BAO's along the line-of-sight (LOS), 
allowing a direct measurement of H (z) has been studied. 
Although the BAO detection along the LOS in the LRG 
catalog of the SDSS DR6 survey used in this reference 
has been called into question [35] , in [33] an enhancement 
along the LOS has been observed in the SDSS DR7 MGS 
which is compatible with a BAO feature. Furthermore, 
[33] argued that a direct measurement of H (z) thanks to 
the BAO peak along the LOS might be possible with a 
photometric survey such that a z /(l + z) = 0.003. 

As an example of what can be achieved with the an- 
gular correlation function, we now study the transversal 
and longitudinal BAO's, using the angular and redshift 
correlation, function obtained in the previous section. 



Transverse correlation function. 



The transverse correlation function is obtained from 



equation ( 16 ) by considering two bins at the same mean 



redshift z and varying their angular separation 8 — 02 — 



e{e,z,z) 



dz 1 W{z 1 ,z)b{z x )G{zi) 

x J dz2W(z2,z)b(z2)G(z 2 ) 
^{4> x ,4> 2 ,r{z x ,z 2 ,9)) , (21) 



where 4>\ and (j) 2 are given in equation ( |13[ ), and S in 
equation (17). As said, numerically it is convenient to 
first evaluate the C™( r )' s ! equation (12 1. Since these 



functions only depend on the cosmology (through the 
power spectrum P(fc, z = 0)), they can be stored and 
used efficiently for a fixed cosmology. 

We expect that RSD increases the amplitude of the 
transversal correlation function with respect to the real 
space result. For example, a spherical distribution of 
galaxies in real space is squashed along the line of sight 
by peculiar velocities. Since galaxies are collapsing to- 
ward the center of the distribution, those between its 
center of and the observer receive an additional red-shift 
from RSD, while those beyond the center are blueshifted. 
When observing the transversal correlation function, this 
leads to increase of the amplitude, since over-dense re- 
gions are enhanced. 

Figure [4] confirms that the increase in clustering am- 
plitude due to RSD along the transverse direction is a 



r(z,z,8) [Mpc/h] 
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FIG. 4: The transversal correlation function. Two redshift 
bins with mean redshift z\ = Z2 = 0.7 are considered. The 
window function is a Gaussian with o z /{l + z) = 0.03. Non- 
linearities are approximated with Halofit. RSD is relevant 
for the shape of the correlation function, and, in particular, 
it shifts the location of the BAO peak. 



relevant effect. The transverse correlation function is 
obtained using a Gaussian radial window function with 
<7,j/(l + z) = 0.03, compatible with the photometric re- 
quirements of, e.g., Euclid [33 05]. As a reference, equa- 
tion Q has been used on the top axis to give the co- 
moving scale corresponding to a certain angle. The BAO 
peak at 9 ~ 3.45° corresponds to about 110 Mpc/h. The 
real space curve is obtained by neglecting all the coeffi- 
cients depending on (3(z) in equation (21). The differ- 



ence in the positions of BAO peak is of order 1% be- 
tween real space and redshift space. In Figure [5] we show 




FIG. 5: Transversal correlation functions. The curves 
correspond, from top to bottom, to mean redshifts z = 
0.3, 0.7, 1, 2, respectively. A Gaussian radial window func- 
tion with cr z /(l + z ) = 0.03 is used. Non-linearities are ap- 
proximately considered with Halofit. 

the transverse correlation function for different redshifts. 
The BAO peak is located at 6 = 7.2°, 3.5°, 2.6°, 1.6° 
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for z = 0.3, 0.7, 1, 2, respectively, which consistently 
corresponds to about 110 Mpc/h in terms of the co- 
moving separation, equation As the redshift z in- 
creases, the angle under which we observe the comov- 
ing BAO scale tbao decreases. This angle is given by 
l(z) = (i +z ) L Da (z) = L / r ( z )> where r(z), given in equa- 
tion ([7]) , monotonically increases with redshift (note that 
this is not the case for Da (z) that shows a turnover at 
z f=3 2). The peak of the curve corresponding to z = 0.3 
is smeared out due to the logarithmic scale, which is used 
because of the large change in the amplitude of the cor- 
relation functions. 



B. Longitudinal correlation function 



d(F 2 )-d(z-!) [Mpc/h] 
50 150 250 
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The longitudinal correlation function is obtained by 
fixing 6 — <f>2 — 4>\ = 0: 



?(0,zi,z 2 ) 



dz-tWiz^z^biz^Giz!) 



x / dz 2 W(z 2 , z 2 )b{z 2 )G(z 2 ) 
xS(0,0,r(z!,z 2 ,0)) , (22) 



where 



2(0,0,r)= ^2 a n±r 



(23) 



711,712=0,1,2 



We expect different effects from RSD on the longitu- 
dinal correlation function than on the transversal one. 
Along the line of sight the over-densities are enhanced, 
but they are also squashed to smaller scales. Instead, 
along the transversal direction, the enhanced over-dense 
regions still extend over the same scales as in real space. 
The squashing in redshift space implies that, along the 
line of sight, the amplitude of the correlation function on 
large (small) scales is reduced (enhanced) with respect to 
the real space result. Figure [6] shows that, for our cos- 
mological parameters, the correlation function becomes 
even entirely negative on scales larger than Az ~ 0.01, 
corresponding to about 20 Mpc/h at z = 0.7. Note that 
on smaller scales non-linear effects becomes important, 
and to compare with observation a treatment of the non- 
linear the 'fingers-of-god' effect is necessary. 

Fig. [6] is obtained from equation ( 22 ) . We set Z\ — 
z — Az/2 and z 2 — z + Az/2 and plot £ s (0, £1,22) as 
function of Az for z = 0.7. A Gaussian radial window 
function with er z /(l + z) = 0.001, compatible with spec- 
troscopic requirements of Euclid [13 H5], is employed. 
The functions (,™{r) are calculated with the specifica- 
tions given in Fig. [4j and the cosmological parameters are 
also the same. As a reference, on the top axis equation ^ 
has been used to show the comoving galaxy separation 
corresponding to a certain redshift separation. The BAO 
peak at Az ~ 0.055 corresponds to about 115 Mpc/h. 
The BAO scale is only very mildly affected by RSD. The 



FIG. 6: The longitudinal correlation function multiplied by 
Az 2 is shown. The redshift Z\ = 0.7 is fixed and we vary 
Az = z%— z\. A Gaussian radial window function with a z /(l+ 
z) — 0.001 is used. Error bars along the abscissas correspond 
to ±<7z. Non-linearities are approximately considered with 
Halofit. 



difference of the peak position with and without RSD 
as shown Fig. [6] is less than 1%. For better visibility 
of the BAO peak we plot Az 2 £ s (Az). It is clear from 
the figure that a redshift resolution which is significantly 
better that a z ~ 0.03(1 + z) is needed to resolve the 
BAO feature. Therefore photometric redshift with the 
accuracy proposed for Euclid would not suffice. How- 
ever, if photometric redshift can be improved to a level 
of <7 2 ~ 0.003(1 + z) as proposed in [44], one might be 
able to see also the longitudinal BAO's with photomet- 
ric redshifts. In Figure [7] we compare the longitudinal 
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FIG. 7: Longitudinal correlation functions with different 
mean redshifts. The curves, with the BAO peak from the left 
to right, correspond to initial redshifts z\ = 0.3, 0.7, 1, 2, 
respectively. A Gaussian radial window function with a z /(l + 
z) = 0.001 is used. Non-linearities are approximately consid- 
ered with Halofit. 



correlation functions at different redshifts. The BAO's 
peak is located at Az = 0.044, 0.054, 0.064, 0.106 for 
Z\ = 0.3, 0.7, 1, 2, respectively, which correspond to 
about 115 Mpc/h in terms of the comoving separation 
r{z2) — r{li). The BAO peak shifts to higher Az = z-i — Z\ 
as the redshift z, at which the transversal correlation 
function is computed, increases according to equation 

0- 

It is also interesting that the longitudinal correlation 

function is fully negative for Az ~ 0.01 and when mul- 
tiplied by (Az) 2 it becomes nearly redshift independent 
before the raise to the acoustic peak. As said before, 
the sign is due to the fact that in longitudinal direction 
clustering is severely squashed due to redshift space dis- 
tortion as it is visible also in direction dependent plots 
which can be found e.g. in Ref. [HJ figure 3] or [JH figure 
1]. On scales larger than the small redshift difference to 
which clustering is squashed, correlations along the line 
of sight are negative. 

The BAO scale inferred from the longitudinal corre- 
lation function is L ~ 115 Mpc/h. This is nearly 5% 
larger than the value expected from the transversal cor- 
relation function, L± ~ 110 Mpc/h. This systematic 
difference is mainly due to the redshift selection func- 
tion. If the redshift bin is relatively wide, the transverse 
BAO scale is reduced by projecting L transversal to the 
line of sight to L± < L. This effect can be considerable. 
L± converges to the longitudinal scale L for a z — > 0, see 
Fig. [8] To take into account that at finite redshift reso- 
lution the BAO peak is not fully transversal, we correct 
the corresponding scale L by adding a small longitudinal 
component L\\ ~ SzH~ 1 (z). 



L 



Sz 



i L {z)D A {l + z) 



Az L /H(z) , (24) 
(25) 



This takes into account an averaged longitudinal contri- 
bution. Solving this expression for F(z) we obtain 



F(z) 



Sz 
~Az~ L 



Az L (z) 

F AP (z) 
F c A I(z). 



1 - 



1 / Sz 



2 \Az L 



2 \Az L 



(26) 



Here we have set Sz = y/jcr z . Naively one might expect 
7 ~ 1, however, we have found that the correction is 
significantly smaller and requires only 7 = 1/8, see top 
panel of Fig. [9] . We suggest that this comes from the fact 
that longitudinally separated galaxies are actually anti- 
correlated on this scale and therefore the contributions 
from galaxies with significant longitudinal separation are 



suppressed in the positive transversal correlation func- 
tion. However, we have not found a satisfactory deriva- 
tion of the pre-factor 7 = 1/8. 

In the top panel of Fig. [9] we compare this redshift- 
corrected Alcock-Paczyhksi function, F AP (z) with 7 = 
1/8 and the uncorrected function F AP (z). Also the real- 
space result is shown, which is even larger. Clearly, the 
difference from the input function F{z) is significantly 
reduced by this redshift-correction. From eq. (24 1 it is 
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FIG. 8: The angle corresponding to the BAO's peak in the 
transverse correlation function is shown for different redshifts 
2 in function of the parameter <r z /(l + z), which determines 
the width of the radial selection function. As a reference, on 
the right it is also indicated the comoving distance according 
to the ACDM model. 

clear, the smaller a z the smaller the difference between 
the longitudinal and the transversal BAO scale. To re- 
move the effect of the redshift bin and check for other 
systematics, we plot in Fig. [9] (lower panel) the real space 
result for F AP (z) / F(z) with redshift bins of zero width, 
i.e. using £ s (6, Z\, 22) without RSD terms to determine 
the BAO peak position. We have also found that it is 
crucial to choose symmetric longitudinal redshift differ- 
ences, £ gal (CJ, z — Az/2, z + Az/2), where z is the redshift 
at which the transverse correlation function is computed, 
and we vary Az. Other choices, e.g. £ gal (0,z, z -I- Az), 
do not lead to a consistent result because the longitu- 
dinal and transversal correlation functions compute the 
BAO's at slightly different redshifts z. Furthermore, we 
have found that the peak position is significantly im- 
proved by using good angular resolution. We compare 
the result of the linear power spectrum and the one using 
halofit to take into account non-linear effects. Finally 
we add formal error bars using a redshift resolution of 
SAz — 0.001(1 + z) from the Euclid spectral survey and 
an angular resolution of 86 = 0.05° and 0.02°. Clearly, 
for z > 1 the angular resolution significantly affects the 
errors. Assuming these formal errors, and allowing for 
several redshift bins, a large survey such as Euclid should 
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FIG. 9: Top panel: z-space (solid) and real-space (long- 
dashed) AP consistency test for cr 2 /(l + z) = 0.03. The 
short-dashed line takes into account the correction F^ t with 
7 = 1/8. Bottom panel: Result without radial window func- 
tion (solid line). Solid error bars consider angular resolution 
80 = 0.02°, and dotted error bars 59 = 0.05°. The dashed 
line consider a linear P(k). For illustrational purposes, the 
analysis is extrapolated up to z tr = 4. 



enable us to measure F{z) with a few percent accuracy. 



IV. CONCLUSIONS 

In this work we have introduced the directly measur- 
able redshift dependent angular correlation function. We 
argue that this function contains the full three dimen- 
sional information about galaxy clustering. The stan- 
dard spatial correlation function requires the assumption 
of a cosmology and a determination of cosmological pa- 
rameters with it can therefore at best be considered as a 
consistency test. We have also found that the shapes of 
the longitudinal and transverse correlation function are 
not only very different from each other, but also very 
different from the real space correlation function. While 
the integral of the latter has to vanish, so that £(r) < 
for large r, we have found that £(0, z,z) > V 9 and 
£(0, z - Az/2, z + Az/2) < V Az > 0.01. 



We used the redshift dependent angular correlation 
function to determine the transverse and the longitudi- 
nal acoustic peak positions as function of the redshift. 
We have proposed to use these functions to perform an 
Alcock-Paczyriski test. We have fully taken into account 
redshift space distortions. Even though these modify sig- 
nificantly the shape of the correlation function, they are 
not very important for the position of the acoustic peak. 

We have found that photometric redshift determina- 
tions with a z /(z + 1) ~ 0.03 are sufficient to deter- 
mine the BAO peak in the transverse correlation func- 
tion, but in order to cleanly determine the position 
of the longitudinal peak, spectroscopic precision with 
a z /(z + 1) — 0.001 are needed. Furthermore, we have 
seen that the peak position of the transverse correla- 
tion function slightly depends on the window function 
for the redshift determination. To achieve an accuracy 
better than about 6% for the transversal BAO scale, one 
has the employ the correction suggested in this work. 
This, together with the correctly symmetrized longitudi- 
nal correlation function allows to obtain an estimation of 
F{z) = (1 + z)H{z)Da{z) which is accurate to about 2%, 
when performing the Alcock-Paczyriski test. 
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Appendix A: Wide-angle RSD using the tripolar 
expansion 

In this appendix we derive expressions for the param- 
and b nino used in the calculation of the cor- 



relation function via Eq. (17) 



We use the coordinates introduced in section |TT] and 
shown in Fig. [T] We can characterize the correlation 
function, equation , by the two directions ni , n 2 from 
the observer and the comoving distance between the two 
galaxies r = ri r 2 = rn. As pointed out in [47l [48] . 
the correlation function in redshift space depends on the 
triangle formed by the observer and the two galaxies. 
This triangle is invariant under rotation about the ob- 
server and hence can be described by one size and two 
shape parameters. For the former we use the separation 
of the galaxies r, while for the shapes we choose the an- 
gles 0i, (f>2, e.g., between ni, n 2 and n, respectively, see 
Fig. [T] We suppose that 4>\ < (f>2 and the distances are 
r*i = r{z\) and r 2 = r(z 2 ). For fixed redshifts Z\ and z 2j 
the distance r is constrained by \r\ — r%\ < r < r% + r 2 . 
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We can evaluate the galaxy correlation function as given 
in Eq. (TTTb, 



£gai(ni,n 2 ,r) 



b 1 G 1 b 2 G 2 



d 3 k 
(2tt) 3 



P(k)e 



k r 



Pi ,2fk 
3 3 

P2 , 2/3 2 



P 2 (n 1 -k) 



PiK-k) 



3 P 2 (n 2 .k) + ^^P l( n 2 .k) 



(Al) 



The most direct way would be to repla ce t he occur- 
rences of k in the integrand of equation (Al) in terms 
of the gradient V r . This is indeed straightforward in the 
plane-parallel limit, because Fourier modes are eigenfunc- 
tions of the plane-parallel distortion operator, that cor- 
responds to the case in which the terms of equation ( 11 ) 



multiplied by Legendre polynomials are negligible, and 
can be written in terms of the inverse Laplacian V~ 2 |47j . 
However, it is non-trivial to obtain an expression for £ ga i 
beyond the plane-parallel limit proceeding in this way. 
In fact, the computation would lead to a large number 
of terms that, although straightforward, are cumbersome 
to evaluate numerically. Despite the non-locality of the 
power spectrum including the wide-angle RSD, the corre- 
sponding correlation function has still been calculated in 
spherical coordinates and expressed in a closed form, e.g., 
in |48j by the introduction of a convenient spherical ten- 
sor and in the approximation of small redshifts z < 1. 
This work has been generalized by [371 SS] to the case 
of arbitrary galaxy redshifts z\, z 2 also for non-flat cos- 
mologies (the a-term there is a selection function) . Even 
with these results, however, the symmetry of the prob- 
lem, namely the dependence of £ ga i on the rotationally in- 
variant triangle determined by {n 1 ,n 2 ,n}, suggests that 
a more natural way to proceed is to expand the expres- 
sion in tripolar spherical harmonics. In particular, we 
use the functions defined in [6l [7J 



%^(ni,n 2) n) 



E 

mi ,777,2 ,777 



t x i 2 i 
mi m 2 m 



x Cg imi (ni)Ce 2m2 (n 2 )Cg m (n) , 



(A2) 



where C lm (n) = 



4t t 
21- 



jYg m are conveniently normal- 



ized spherical harmonics and ( „ x ^ 2 2 m) ls the Wigner 
3-j symbol. These function form an complete orthogo- 
nal basis for expanding spherically symmetric functions 
depending on three vectors, which are invariant under 
global rotations. In [SJ [7J it has been shown that, in the 
approximation z <C 1, this expansion leads to a more 
compact form for £ ga i than those obtained in previous 
works. We generalize the work of [SI [7] by allowing arbi- 
trary galaxy redshifts z\, z 2 and a generic function a(z). 
After some algebra one finds that the only non- 



vanishing coefficients of the expansion 
£ g ai(rii,n 2 ,n,r) = b 1 G 1 b 2 G 2 

x P W (^i^2)S w ( ni ,n 2 ,n) , (A3) 



are 



B 220 (r,z 1 ,z 2 ) = -^=(3 1 p 2 ( 2 (r), 



B 222 (r,z 1 ,z 2 ) = ^^ 2 Cl(r) 
B™(r,z 1 ,z 2 ) = ^§p 1 p 2 ti(r), 



35 



B™\r,z u z 2 ) 



B 022 (r,z u z 2 ) = - Q/3 2 + -ft/fc) V5 C 2 2 (r) ,(A4) 

which are independent of the angles (as they do not in- 
volve the a-terms), and 

B^r,z 1 ,z 2 ) = - a ^[ Pl + ^ j tf(r), 
B^r,z 1 ,z 2 )= ( ^*(p 2 + ^) Cl(r), 



B 121 (r, Zl ,z 2 ) 
B 123 (r, Zll z 2 ) 



r 2 V 3 



n \/T5 



B 2n (r, Zl ,z 2 ) = - 
B 213 (r, z\, z 2 ) = - 



n V15 

a 2 U 2 2V2 
r 2 

a 2 n 2 2V7 



r 2 V15 



Aft Ci(0 , 



B™(r,zi,Z2) = ~ aiHia2H2 ^ 2 Co °(r) , 



B 112 (r )Zl ,z 2 ) 



T\T 2 ^3 

OL X U x a 2 U 2 /10 o 

V T^i & C 2 ( r) 

r\r 2 V 3 



(A5) 



If one consider f3 x = fi 2 = /, i.e., neglecting the bias and 
setting z\ = z 2 , and a% = a 2 = 2/T-L, this result agrees 
with [7J . We also verified the consistency with (TO] . 

To proceed further we must choose a system of coordi- 
nates. As shown in [5], however, once the coefficients of 
the tripolar expansion have been calculated, it is trivial 
to write the correlation function for a given coordinate 
system. In fact, this step involves only simple algebra 
that, although tedious, can easily be carried out using a 
computer algebra package. 

We use the coordinate system a) discussed in [B] since 
it gives the most compact form for £ ga i. As shown in 
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Figure [Tj the direction z is orthogonal to the plane that 
contain the triangle formed by the observer O and the two 
galaxies. Hence, all the vectors r l5 r 2 , r have latitude $ = 
7r/2. The galaxy separation vector r has zero longitude 
= 0, i.e., it is parallel to the abscissas. 

With this choice of coordinates we can evaluate the 



elements of the tripolar basis, equation ( A2 ) , which so far 
have been written in terms of the re-norma ized spherical 
harmonics CV m (n). Using the notation n = {$,0}, we 
have 

S^^ni, n 2 , n) = S tl e 2 t({n/2, x } , {tt/2, 2 } , {tt/2, 0}) • 

Evaluating these expressions, the summation over £i,£ 2 
and I in Eq. ( A3 1 , using ( A4 1 and ( A5 1 for the non- 



vanishing functions B 1 2 leads to 

£gai(si, ^l) z 2 ,4>2,r) = b(z 1 )G(z 1 )b(z 2 )G(z 2 ) 
> </„,„, cos (n.i0i) cos (n 2 0a) 



ni ,ri2— 0,1,2 



+6 ni „ 2 sin(n 1 1 )sin(n 2 2 )] . 

(A6) 

The coefficients a ni „ 2 and 6 ril , i2 can be easily calculated 
with a numerical algebra package. We first write all the 
non-vanishing coefficients which do not involve a-terms: 



l + J(/3i+/3 2 ) + ^/3i/3 2 )^(r) 



2 

21' 



£(A+/9 2 ) + ^i^)C 2 2 (r) 



+ 140 /?1/32 C|W 



O20 

«02 

022 
l>22 



- ( ^2 + ^/3i/3 2 ) C!(r) + ^ftft Cl(r) , 



1 



I 19 



^i/3 2 Co 2 M ' ' 2< ■ 1 ' ' J 



(A7) 



The coefficients that involve a-terms are: 



aio 

aoi 
an 

021 

ai2 
6n 

&21 

bi2 



aiKxPx + ^atflM^) -Ci(r) 

5 / n 

ai^i /3i/3 2 j . , 

"C 3 M » 
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OL 2 U 2 p 2 + \a. 2 U 2 p x p 2 \ 

5 / r 2 



. «2%2 /3i/3 2 u 

aihLxa 2 M. 2 PiP 2 - 

o so l r / 

3 r x r 2 



2u{H.\a 2 'H.2 PiPi a0 



'"2 



rir 2 



C 2 » 



a 2 H 2 PifafX,^ , 3 o; 2 , H 2 P\P 2/hl 
10 r 2 



C 3 x (r) 



-CiW 



-C 3 ( r ) 
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ai'Hia 2 , H 2 /3i/? 2 , 
3 rir 2 
a 1 T-L 1 a 2 U 2 PxPi a0 



rir 2 



C 2 >) , 



a 2 H 2 I3xl3 2 x a 2 H 2 P\P 2 x , , 
-Ci(r) — C 3 U) 



'2 



'2 



-Ci M 



''1 



ai^i /3i/3 2 x 



(A8) 



Setting Pi = P 2 = f and ai = a 2 = 2/H, we re- 
cover the results of [Tj- We recall that we are assuming 
01 < 02 (Figure [I]), this constraint can be inverted us- 
ing the symmetry £gal(n 2 , m, n, r) = £ ga i(-ni, -n 2 , n, r) 
evident from equation (111. 
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